A QCD Sum Rule Study of 9+ in Nuclear Matter 
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We consider a [ud\^s current, in the finite-density QCD sum rule approach, to investigate the 
scalar and vector self-energies of the recently observed pentaquark state O"'"(1540), propagating in 
nuclear matter. We find that, opposite to what was obtained for the nucleon, the vector self-energy 
is negative, and the scalar self-energy is positive. There is a substantial cancellation between them 
resulting in an attractive net self-energy of the same order as in the nucleon case. 
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The possible existence of a narrow exotic baryon 0+ 
with strangeness -t-1 is, nowadays, possibly one of the 
most exciting topics in nuclear physics. There is a large 
number of experiments with evidence for the existence 
of this state [l| , and a similar number of high-energy ex- 
periments that see no evidence 5^. Even lattice QCD 
simulations of pentaquarks by several groups have not 
converged yet 0|. The apparent contradiction be- 
tween evidences for and against the existence of 8+ 
might be resolved if there is a production mechanism 
which is present in some experiments and absent in oth- 
ers. One such production mechanism was proposed in 
refs. 0,01 and is related with the formation of the cryp- 
toexotic A^*(2400) resonance, which would decay into 
6+ K-. The TV* (2400) could be formed in the reac- 
tions 7 + p ^ 7r+ + iV*(2400) and n- + p ^ N*{2A00) 
p|, which arc not present in e^e^ collisions and could be 
dificult to happen at high-energy collisions. Other pos- 
sible production mechanism is meson exchange between 
the nucleon and another hadron. Since at high energies 
all exchanges with the non-vaccum quantum numbers die 
out and only the gluonic pomerons survive, it may be dif- 
ficult to excite the pentaquark by soft gluons on nucleon. 



If 0+ is indeed produced through the decay of the 
A^*(2400) resonance or through meson exchange, then 
one should expect that it could be seen in heavy ion colli- 
sions at RHIC, since these two processes could happen in 
the collision of nuclcons and many "comoving" hadrons 
produced in the collisions at RHIC. However, if this is 
the case, the O"*" would be formed in a nuclear medium 
which could change its mass and decay width. There- 
fore, to be able to identify positively the 0+ signal in 
heavy ion collision at RHIC, it would be very important 
to know how the nuclear medium affects the O"*" char- 
acteristics. In this work we will use the QCD sum rule 
(QCDSR) approach in nuclear matter [1,0, to study 
the pentaquark 0+ in a finite density medium. 

The finite density QCDSR approach focuses on a cor- 
relation function evaluated in the ground state of nu- 
clear matter, instead of the QCD vaccum (as in the usual 
sum rules). For spin-i baryons, this function can be de- 
composed into three invariant functions of two kinematic 



invariants. The appearence of an additional invariant 
function, compared with the vaccum case, is due to an 
additional four-vector: the four-velocity of the nuclear 
medium, u^, which, together with the the nuclear den- 
sity, pn, characterizes the nuclear matter ground state. 
The quasibaryon excitations are characterized by scalar 
and vector self-energies. By introducing a simple ansatz 
for the spectral densities, one obtains a phenomenological 
representation of the correlation function. 

The correlation function can be also evaluated at large 
spacelike momenta using an operator product expansion 
(OPE). This expansion is written in terms of the matrix 
elements of composite quark and gluon operators, evalu- 
ated in the nuclear matter ground state: the in-medium 
condensates. By equating the OPE and phenomenolog- 
ical representations of the correlation function, one ob- 
tains QCD sum rules that relate the baryon self-energies 
to the in-medium condensates. 

Several zero-density QCDSR investigations of 8+ al- 
ready exist where the authors have used different inter- 
polating fields 0, 0, 0, Q, Here we use the inter- 
polating field suggested in refs. [3, 0| : 

Vi^) = eabcedef^cfg[Ua{x)Cdb{x)][u^{x)Cj5de{x)]Cs'^{x), 

(1) 

where a, b, c, ... are color index and C — — is the 
charge conjugation operator. In Eq. each diquark pair 
has spin and isospin zero and is in the 3 representation of 
color SU(3). The total current has isospin zero, positive 
parity and spin 1/2. In ref. 13] it was shown that the 
ground state described by this current has negative parity 
and a mass compatible with the experimental 8"'' mass. 

The QCDSR for 8+ at finite density is based on the 
correlation function defined by 



n{q)=t / d4^e*«-^(^'o|T7?(x)^(0)|^o) 



(2) 



where j^E'o) represents the nuclear matter ground state. 

The correlator in Eq. Q can be decomposed in three 
distinct structures ^ 

Uiq)=Usiq^,q-u) + Uq{q^,q-u)^ + Uu{q^,q-u)i. (3) 

In vaccum Hg and H^ become function of q^ only and H^ 
vanishes. For simplicity we will work in the rest frame of 



2 



nuclear matter, which imphes that = (l7 0)- There- 
fore, Ui{q'^,q ■ u) lli{qo, \q\) {i = {s, q, u}). 

In the phenomenological side, the analytic properties 
of n(q) can be studied through a Lehman representa- 
tion, which leads to a dispersion relation in qo, for each 
invariant function, of the form Isl 



n^(qo, \q\) 



1 



- ,^An.(.,|,1) 



(4) 



where we have omitted polynomials arising from the 
contour at large \qo\, which will be eliminated by the 
Borel transformation. The discontinuity, defined by 
AIli{u}, \q\), contains the spectral information on the 
quasiparticle, quasihole, and higher-energy states. 

In vaccum, the spectral weights for baryon and an- 
tibaryon are related by charge conjugation symmetry and 
one usually parametrizes the spectral density as a sin- 
gle sharp pole, representing the lowest resonance, plus a 
smooth continuum, representing higher-mass states. At 
finite density, the ground state is no longer invariant un- 
der charge conjugation and, therefore, the spectral den- 
sities for baryon and antibaryon are not simply related. 

The width of 8+ in free space is very small and can 
be ignored on hadronic scales At finite density, the 
width of 9+ can be broadened due to strong interac- 
tions. Since the introduction of one extra parameter in 
the spectral density would reduce the predictive power of 
the sum rule, here we assume that a sharp pole hypotesis 
is still valid at finite density. In the context of relativis- 
tic phenomenology, we assume that 0+ couples to the 
same scalar and vector fields as the nucleon and the hy- 
perons in nuclear matter. Therefore, we parametize the 
discontinuites as ISl 



An,(c^,|<j1) = 

AIlg{u;,\q]) = 



An„(^,|Ql) 



where 
and E„ 



2E* 



[Siio - E,) - Siio - E,)] , (5) 



me 



E* = 



q^, Eq = Y.^ 
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jq — - . /\Q measures the coupling of the 
interpolating field with the physical in the medium 
and has opposite signs in Alls depending on the parity 
of 0+ 17]. The scalar and vector self-energies of 0+ in 
nuclear matter are given by and respectively. The 
positive- and negative-energy poles are at Eq and Eg re- 
spectively. In Eq. (jsj we have omitted the contributions 
from higher-energy states, which will be included later 
in the usual form, i.e., being approximated by the OPE 
spectral density starting at an effective threshold. 

At finite density, the OPE for the invariant functions 
takes the general form 8, 9] 



n.(go,|'7l) 



^C;(go,k1)(0„)p„, 



(6) 



where (On)p„ — {^a\On\'^o), are the in-medium conden- 
sates. The Wilson coefficients, C^(go, M) depend only on 
QCD Lagrangian parameters. Therefore, all the density 
dependence of the correlator is included in the in-medium 
condensates. 

Following ref. we separate the invariant functions 
into two pieces that are even and odd in qo: 



n,(go, m = Ufiql m + qoTlfiql, \q]) , 



(7) 



In ref. [Sj it was shown that the choice go — Eq completely 
suppresses sharp excitations at Eq and also strongly sup- 
presses a broad excitation in this vicinity. Since we are 
interested in the positive energy pole, we will use og = Eg. 

By using the quark propagators given in refs. 
and working to leading order in perturbation theory, up 
to dimension 5, we get (after Borel transforming both 
sides of the sum rules) : 
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EqM'^E^ f^^^. 

237r64!5! ^'^ ' 2%63!5! 
1 , . , A 4(5ssV-gs)_^^^g^ Eq{M^E2 + 2<tM^E^) 



7r64!6! 



267r63!5!3 
Eq{gsq'^(J-Qq) 



227r64!5! 
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237r64!5! 



M^E2 , 



(8) 
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^((^^^^o.),„-i(gt^Z?o<z),„) , (10) 
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where we have defined 

fe=0 

which accounts for the continuum contribution with sq 
being the continuum threshold. 

To extract the self-energies from the above finite den- 
sity sum rules, one has to know the values of the in- 
medium condensates. To first order in the nucleon den- 
sity, one can write {O) ~ (O) -|- {0)nPn where {0)n 
is the spin-averaged nucleon matrix element. The sim- 
plest in-medium condensates are {q''q)p^ and (s^s)p„- 
Since the baryon current is conserved, {q^q)p^ and 
{s^s)p^ are proportional to the nucleon and strangeness 

densities: {q^q)p^ — ^Pn and (s^s)p„ = 0. These 
are exact results. We extract the values of the other 
condensates from refs. d, 0, 0|: (gg)p„ = (gg) + 

{g^,G^)~47:^{0.65GeY)pN, {q^iD^q)p, = (0.18GeV)pw, 
{sHD,s)p^-^{ss)p^ = (0.018GeV)p^, {9sqcr-Qq)p^ - 
^l{m)p,,.{9sS(J-gs)p^ = ml{ss)p^, {qiDoiDoq)p^ + 
U9sqa-gq)p^ = i0.3GeV^)pN, (siDoiDos) + 
l{9sS(J-Qs)p^ = y(0.3GeV2)/3Ar, {gsq^o-Gq) p^ = 
(-0.33GeV^)pAr, {gss'^a-gs)p^ = y{9sq^ (^-Gq) p^, 
(qHDoiDoq) p^ + ^{gM^a-gq)p^ = (0.031 GeV^)^^, 

{sHDoiDos)p^ + ±{g,s^a-gs)p^ = y{OmiGeV^)pN 
where is the nucleon a term and y — (ss) at / (gg) at is 
a parameter that measures the strangeness content of the 
nucleon. In this work we use cr at = 45 MeV, rriu + rud = 
12MeV and y in the range < ?/ < 0.5 [l^- The values 
used for the strange quark mass and vacuum conden- 
sates are the same as used in ref. nis = 0.13 GeV, 
(qq) = -(0.23)3 GeV^ (ss) = 0.8{qq), ml = 0.8 GeV^ 
and (5sG^) = 477^(0.33 GeV)^. Nuclear matter satura- 
tion density is taken to be pn = Po = (110 MeV)'^. For 



the value of the three- momentum we follow refs.|8ll9lllCl| 

and fix it at |g| = 270 MeV. 

We evaluate our sum rules in the range 1.7 < < 
2.5 GeV^ As in ref. ^ we found that the left-hand side 
of Eq.@ is negative indicating that the parity of 0+ is 
negative [l7| . Since the sum rule in Eq. (|5J is very unsta- 
ble, in this work we neg lect it. To eliminate the coupling 
Aq^, we follow refs. |llLll2| and derive another sum rule 
by taking the derivative of Eq.© with respect to M"^ 
and dividing it by Eq. ©. We also divide Ea. l(Tn|) by 
Eq.® obtaining two sum rules that are independent of 
Aq^. To quantify the fit of the left- and right-hand sides 
of the sum rules and to determine the self-energies, we 
apply the logarithmic measure defined in refs. aa,UJ- 
To check the validity of this procedure we start by analyz- 
ing the sum rule in vacuum. The optimized values for the 
0+ mass and continuum threshold are me = 1.64 GeV 
and So = 3.55 GeV^, which are in a good agreement with 
the values found in ref. [isj . 

In Fig. 1, using y = 0.3 |l8| |. we show, in the upper part, 
the left- (dotted line) and right-hand sides (dashed and 
dash-dotted lines) of the sum rules (in units of GeV^) as 
functions of . We see that there is a good agreement 
between them. In the lower part of Fig. 1 (solid line) we 
show the effective 0+ mass. The optimized values for the 
effective mass, the vector self-energy and continuum 
threshold are: uiq = 1.75 GeV, E„ = -150 MeV and 
So — 3.6 GeV^. The first astonishing result is that, op- 
posite to what was obtained in the nucleon and hyperon 
cases, the scalar self-energy is positive and the vector 
self-energy is negative. Therefore the effective 0+ mass 
in medium, ttiq = ttiq + Sg, is bigger than the free Q'^ 
mass. However, since Ss — 110 MeV, there is still a 
substantial cancelation between Es and E„ in medium. 
Looking at Eq. (fTn|l it is easy to understand why we get 
< 0, since, in contrast to the nucleon and hyperon 
cases, all the terms in the right-hand side are negative. 
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FIG. 1: The left- and right-sides of the sum rules (upper lines) 
and the effective mass (solid line) as functions of the Borel 
parameter M'^, for y — 0.3. 



The origin of this difference in the sign can be attributed 
to the existence of an antiquark in the interpolating 
field. Since and S^, are essentially the real parts of 
the optical potential, this qualitative result shows that 
there is an overall attractive 0— nucleon interaction. 
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suits are very sensitive to the value of y. In the range 
< y < 0.5 we got 50MeV < < 150 MeV and 
-90MeV > > -190 MeV. What is very interest- 
ing is that the sum Eg -I- S^,, which can be associated 
with the depth of the O"*" potential in nuclear matter, 
remains independent of y and is about —40 MeV. How- 
ever, this value is very sensitive to the value of the gluon 
condensate. If we change the value of the gluon conden- 
sate to {g'^G'^) = 0.24 GeV^, as used in ref. we get 

(for y = 0.3): E^ - 90 MeV and E.„ 180 MeV which 

would imply in a potential depth of about —90 MeV. 

The depth of the potential in nuclear matter, U, 
was studied in ref. |l9l | using a relativistic mean field 
framework, where couples with scalar and isoscalar- 
vcctor mesons. They found -90 MeV < U < -45 MeV 
depending on the values used for ■m*^/mN (the ratio of 
the nucleon mass), and the coupling constants. In our 
work we found that the potencial depth depends strongly 
on the value of the gluon condensate. In the range 
0.47 GeV^ < {g^G^) < 0.24 GeV* we got -40 MeV > 
U > — 90 MeV, which is compatible with the findings of 
refTli. 

It is also interesting that, in spite of rriQ being bigger 
than rriQ, the energy of the quasi-0 in nuclear matter, 
Eq, is smaller than me and it is the most stable result 
of our calculation. It does not depend on y neither on 
(ggG^). While me, ttiq and E.y varies significantly with 
y and {g'lG'^) we got Eq ~ 1.60 GeV for aU the values of 
the parameters (the smallest value obtained for the 0"*" 
mass was me — 1.64 GeV). 
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FIG. 2: The ratio m@/me (solid line), the vector and scalar 
self-energies (dotted and dashed lines respectively) as func- 
tions of the Borel parameter M^, for y — 0.3. 

The optimized results for the ratio rriQ/mQ (solid 
line) and for the scalar (dashed line) and vector (dot- 
ted line) self-energies as functions of are plotted in 
Fig. 2 for y = 0.3. We see that the curves are quite 
flat, indicating a weak dependence of the predicted re- 
sults on M^. The results are also not sensitive to the 
value of \q\. Using |gl = instead of \q\ = 270 MeV 
does not alter our results. On the other hand, the re- 



-0.5 ' • ' ■ ' ' ' • ' ' ' • ' • ' ^ 

0.2 0.4 0.6 0.8 1 1.2 1.4 

P/Po 

FIG. 3: The ratio rriQ/me (solid line) and the vector self- 
energy (dashed line) as functions of the baryon density, for 
y = 0.3. 

All the results given above were obtained at the nuclear 
matter saturation density po. In Fig. 3 we show the den- 
sity dependence of the scalar and vector self-energies. We 
see that most of the variation happens for < 0.5po- 
Since we are expanding the density dependence of the 
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condensates up to first order in the nuclear density, the 
obtained density dependence of our results is compatible 
with the approximations made. 

We would like to point out that, since our sum rules 
do not depend on the four-quark condensates, our results 
are free from the uncertainties associated with the density 
dependence of these condensates, which were the biggest 
source of uncertainty in the case of the nucleon and hy- 
perons studied in refs. |1,0>0|- In that case, the scalar 
self-energy was very sensitive to the density dependence 
of the four-quark condensates, and one could even get 
a repulsive net self-energy. The results were a bit more 
stable for the vector self-energy. Comparing the values 
of the vector self-energies one has from refs. 0, 0, llfij |: 
E^/Ef ~ 0.3 - 0.4, E^/E^ ~ 0.8 - 1.1. And here we 
got |E®|/E^ ~ 0.3 — 0.5. In terms of relativistic hadronic 
model these results would imply that the coupling of the 
hyperon E to the Lorentz vector field is very similar to 
the corresponding nucleon coupling, while the coupling 
of the hyperon A is similar to the coupling of 8 and both 
are much weaker than the correponding nucleon coupling. 
These results can be understood in terms of the interpo- 
lating fields used to study these baryons. For A there is a 
[ud] diquark with spin and isospin zero plus the strange 
quark, that carries the spin of A. Assuming no admix- 
ture of strange quark content in the vector meson lo (the 
Lorentz vector field), there will be no coupling between A 
and Lu. The same happens in the case of (see Eq. (Q). 
The situation is different in the case of E, since in its in- 
terpolating field there is a [qs] diquark with spin zero. 
Therefore, it is the light quark that carries the spin of E, 
as in the nucleon case. 

As a final remark, we would like to mention that in 
ref. 1^ it was suggested that if the parity of 0+ were pos- 
itive, then there would be a strong attractive 0— nucleon 
interaction. In our calculation we got a negative parity 
for 0"*". Its interaction with the nucleon is still attractive, 
and roughly of the same order of the nucleon-nucleon in- 
teraction. Bigger attraction is obtaned with smaller val- 
ues of the gluon condensate. 
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